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Abstract 

We present a systematic study of dijet suppression at RHIC using the VNI/BMS parton cascade code. We examine 
the medium modification of the dijet asymmetry Aj and the energy distribution within the dijets. Understanding the 
sensitivity of these observables to properties of the deconfined medium and to experimental factors is vital if dijets 
are to be useful for QCD tomography. VNI/BMS provides a controllable test-bed with sufficient complexity to model 
jet modification without confounding results with fluctuations from hydrodynamics and hadronization. Dijets are 
examined under the modification of: the jet transport coefficient q; the path length of leading and sub-leading jets; jet 
cone angle and the jet-medium interaction mechanism. We find that Aj is very sensitive to the distance traveled by the 
secondary jet, the jet shape is dominated by q and the nature of the interaction mechanism. 

Keywords: parton cascade, jet suppression, quark-gluon plasma, dijets, Hard Probes 2012 



1. Introduction 

The suppression of high energy E, ~ 100 - 200 GeV dijets in heavy ion collisions has recently been observed at 
the LHC (HID. These results have shown the feasibility of using dijets as correlated probes of jet modification in hot 
QCD matter. Proposed upgrades to PHENIX along with the continuing development of jet reconstruction at STAR 
will allow access to dijets in the RHIC kinematic window E t ~ 15 - 65 GeV. The dijet asymmetry is defined as [1] 

Ei f — E f 

&t,t + E t>s 

where E,j is the transverse energy of the leading jet and E ttS is that of the sub-leading jet. The asymmetry distribution 
measured in Pb+Pb collisions at the LHC has been successfully reproduced with a variety of models including the 
one considered here lU|5]|6l|3. Recent work by Renk [18] indicates that these results may be relatively insensitive to 
the fine details of jet energy loss, suggesting that Ay may be unsuitable for tomographic purposes. If this is indeed the 
case, are there other more differential dijet observables for jet measurements at RHIC? 

With these issues in mind we have used the VNI/BMS parton cascade to systematically explore dijet suppression 
at RHIC energy scales under controlled medium conditions. We examine the modification of dijets under variation of 
the medium radius and temperature, the strong coupling constant and the jet definition in terms of the Anti-Kt cone 
angle. The VNI/BMS parton cascade model code J9] [TOl provides access to the full jet/medium development at a 
fixed q. We run the code in a static uniform-medium mode. The medium is modeled as a torus of a given radius. 
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The dijet propagation lengths are generated as chords centered on uniformly sampled hard-collision vertices within 
this torus. The model includes a partonic medium which is treated on an equal footing with the jet. This allows jet 
partons to escape into the medium and vice-verso. All results are presented at the partonic level only. We use anti-kt 
jet reconstruction [ 1 1 1 throughout to give a somewhat realistic treatment of the jet measurement process. 

2. The Parton Cascade 

The parton cascade model (PCM) is a Monte-Carlo implementation of the relativistic Boltzmann transport of 
quarks and gluons 

f Q^Fkix, p) = J] CiF k (x, p). (2) 

i 

The collision term C,- includes all possible 2 — > 2 interactions and final-state radiation 1 — > n 

dF k {x, p) = J" Y] ^;IMf 5 4 (A,, - Pout) D(F k (x, p*)), (3) 

dTj is the Lorentz invariant phase space for the process j, D is the collision flux factor and 5 , is a process depen- 
dent normalization factor. A geometric interpretation of the total cross-section is used to select pairs of partons for 
interaction. Between collisions, the partons propagate along straight line trajectories. 

In the VNI/BMS implementation of the PCM outgoing off-shell partons are brought back on-shell through a 
medium modified time-like branching. The partons created in this process are subject to a Monte-Carlo LPM effect 
fT2"l l6l. The strong coupling constant for scatterings is held fixed at a s = 0.3 although we shall explore the effects of 
its variation below. A quark-gluon plasma is simulated as a box of thermal quarks and gluons generated at some fixed 
temperature. Periodic boundary conditions are imposed on the box, whose size is selected to be large enough that if a 
simulated jet wraps around it will not interact with its own wake. 

The partonic contents of jets, created by the event generator PYTHIA-8 lfT3lfT4 l and evolved down to Qo = 1 GeV, 
are injected into the medium. Each sub-jet is evolved separately. The jet evolution is recorded without reference to 
a particular jet definition (cone-angle, energy cuts, jet reconstruction algorithm, etc). Each particle in the initially 
inserted jet is marked as being "jetty", this jettiness tag is then iteratively applied to all partons that interact with 
an already jetty-labeled parton. The entire time-evolution of the jet can then be reconstructed using a suitable jet 
finder (FastJet [ 15 1) and jet definition. In this way the influence of varying jet definitions can be readily studied. This 
process factorizes the initial vacuum radiation process from the medium induced radiation, this is a weakness that will 
be addressed in future work. 

VNI/BMS is a simple enough jet-suppression model. While it may not be suitable for use as a full event generator, 
lacking hydrodynamical flow and hadronization, the active medium and medium-modified radiation make it a useful 
test bed for many but not all jet modification features. 

3. Results 

We generate candidate Monte-Carlo dijet events from pp collisions at yfs — 200 GeV, anti-kt reconstruction was 
used to identify satisfactory dijet events IfTTII . The following kinematic requirements were always imposed: we fix 
the minimum momentum of an identified jet to be p, jm ; n = 5 GeV, A^nl < 1.1, A0i2 > n/2 where Arjn and A0i2 are 
measured between the two reconstructed jets in the event. The constituent partons of each identified jet are translated 
from PYTHIA into the parton cascade medium for modification. We consider the dijet asymmetry Ay and the partonic 
jet profile. 

In Fig. [T] we show the variation of Aj with the medium radius. The vacuum jet distribution falls off very quickly, 
jets with an initial E t ,t > 65 GeV are extremely rare. Increasing the medium temperature from 250 MeV to 350 MeV 
(see Fig.|2jl leads to a depletion of jets with a small modification, shifting the dijet distribution towards higher Aj 
values. By applying an energy cut to differentiate the leading and sub-leading jets in a dijet pair we have biased our 
results. We select for leading jets which travel a short distance, and are less modified, along with sub-leading jets that 
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Figure 1 : The dijet asymmetry as a function of the path length of 
the leading and sub-leading jets, for a medium with T = 350 MeV. 



Figure 2: The influence of the medium temperature on the asym- 
metry of dijets with E,\ > 20 GeV. 



travel a very long distance and so experience a larger integrated q. This surface bias gives a large contribution to Aj 
relative to the vacuum result. For the remainder of the analysis we have fixed R me d — 5 fm . 

Since q cc r 3 the dijet asymmetry depends strongly upon the medium temperature. In Fig. [2] we show results 
of varying the medium temperature. Increasing medium temperature and therefore q leads to increased jet-medium 
interactions and a strong swing in the observed Aj. 
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Figure 3: The variation of Aj with R the anti-kt jet cone radius. 
The medium radius is fixed at R me j = 5 fmwith T = 0.35 GeV 
and E t i > 20 GeV. 



Figure 4: Variation of the jet radial profile at T = 350 MeV. Lead- 
ing jet profiles are shown as solid lines and sub-leading jets are 
dashed. Results for elastic energy loss only are shown in red; those 
including radiative and elastic energy loss are in orange. 



In Fig. [3] we show Aj as a function of the anti-kt cone angle R for jets in a medium with T = 350 MeV. As R is 
increased the amount of dijet modification is significantly reduced. It is important to note that the only partons which 
can be included as part of the measured jets were either directly created by the jet or are those which have scattered 
with jet partons. The thermal medium is currently artificially excluded from the jet finder, this removes uncertainty 
associated with background removal. As R increases more of the relatively soft radiated partons and forward scattered 
medium partons are included in the jet definition along with original hard core. This leads to higher reconstructed 
jet energies at larger R's, which in turn leads to the observed reduced Aj. The effect is proportional to the distance 
traveled by the sub-leading jet, jets which have traveled shorter distances have built up less of a cloud of soft partons 
and picked up fewer medium partons by elastic forward scattering. The rate of transverse diffusion of these soft 
partons is proportional to q, so the jet-cone/cloud will be wider at higher medium temperatures. 
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Let us now examine the modification of the radial jet profile under the same set of factors. We define the radial jet 
profile as the ratio of jet energy reconstructed within a certain jet-cone radius R relative to the reconstructed energy 
at R = I. This gives a normalized radial profile which is very sensitive to variations in the medium temperature, the 
jet interaction mechanism and the strong coupling. This observable could be experimentally obtained by a similar 
process of iterative jet reconstruction. In Fig.|4]we show the radial profile for both elastic and radiatively modified 
jets. The leading jet profiles (solid lines) are somewhat modified compared to the vacuum jets (black profiles), we 
see very little separation between the radiative and elastic only leading jet profiles. The sub-leading jet profiles are 
dramatically modified compared to the vacuum and leading jet profiles. The elastic and radiative profiles clearly 
separate, the radiative sub-leading jets become broader and softer than the elastic only. Both sets of sub-leading jets 
become much broader and softer compared to the leading jets. 

4. Conclusions 

We have examined the modification of RHIC scale dijets by a pQCD medium over a range of temperatures. We 
find that the dijet asymmetry is quite sensitive to the amount of medium encountered by the jets, and very sensitive to 
variations in q. The leading jets show a strong surface bias and lose very little energy making their E t a reasonable 
proxy for their initial energy. As such A j appears to give a reasonable measure of the modification of sub-leading jets 
at these scales. 

The jet profiles we present show the strong modification of the sub-leading jets as a function of temperature and 
strong coupling. They also differentiate between interaction mechanisms, understanding the nature of the jet-medium 
interaction as well as its strength is of paramount import. Although these profiles may well be sensitive to the non- 
trivial hadronization process, we believe that they could provide a very informative window into the real nature of jet 
interactions with the QGP. 

Dijets at RHIC scales are likely to be strongly modified by the presence of the deconfined QGP medium. The 
observables we have discussed are sensitive to many aspects of this modification and suggest that further jet mea- 
surements at RHIC will provide valuable insights into the nature of the QGP and into the applicability of pQCD jet 
suppression models . 
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